Three fermions in a box at the unitary hmit: universahty in a lattice model 
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We consider three fermions with two spin components interacting on a lattice model with an 
infinite scattering length. Low lying eigenenergies in a cubic box with periodic boundary condi- 
tions, and for a zero total momentum, are calculated numerically for decreasing values of the lattice 
period. The results are compared to the predictions of the zero range Bethe-Peierls model in con- 
tinuous space, where the interaction is replaced by contact conditions. The numerical computation, 
combined with analytical arguments, shows the absence of negative energy solution, and a rapid 
convergence of the lattice model towards the Bethe-Peierls model for a vanishing lattice period. 
This establishes for this system the universality of the zero interaction range limit. 

PACS numbers: 03. 75. Ss, 05. 30. Fk, 21. 45. -|-v 



Recent experimental progress has allowed to prepare 
a two-component Fermi atomic gas in the BEC-BCS 
crossover regime and to study in the lab many of its 
physical properties, such as the equation of state of the 
gas and other thermodynamic properties, the fraction of 
condensed particles, the gap in the excitation spectrum 
corresponding to the breaking of a pair, the superfiuid 
properties and the formation of a vortex lattice, the effect 
of a population imbalance in the two spin-components 
and the corresponding possible quantum phases, . . . 

7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. 

The key to this impressive sequence of experimental re- 
sults is the use of Feshbach resonances [I^l : an external 
magnetic field {B) permits to tune the two-body s-wave 
scattering length a almost at will, to positive or nega- 
tive values, so that one can e.g. adiabatically transform 
a weakly attractive Fermi gas into a Bose condensate of 
molecules. Interestingly, close to the resonance, the scat- 
tering length diverges as a oc — l/(i? — Bq) so that the in- 
finite scattering length regime (|a| = oo) can be achieved. 
When the typical relative momentum k of the particles 
further satisfies kb <^ l,k\re,\ <C 1, where h is the range 
and re the effective range of the interaction potential, the 
s-wave scattering amplitude between two particles takes 
the maximal modulus value fk = —l/{ik): This is the 
so-called unitary regime, where the gas is strongly, and 
presumably maximally, interacting. 

The unitary regime is achieved in present experiments 
for broad Feshbach resonances, that is for resonances 
where the effective range r^ is of the order of the Van der 
Waals range of the interatomic forces [l^, [l^ . Examples 
of s-wave broad resonances are given for ^Li atoms by the 
one at Bq ~ 830 G P, i, d, 0] or also for ^Or atoms at 
Bo ~ 200 G [|. On a more theoretical point of view, the 
unitary regime has the striking property of being univer- 
sal: E.g., the zero temperature equation of state involves 
only h, the atomic mass m, the atomic density and a nu- 
merical constant independent of the atomic species; this 
was checked experimentally, this also appears in fixed 



node Monte Carlo simulations [20|, |2l| and more rece ntly 
in exact Quantum Monte Carlo calculations [12, [l^, |24| . 

In Refs. [l^, I exact Quantum Monte Carlo simula- 
tions at the unitary regime are performed using a Hub- 
bard model. From the condensed matter physics point of 
view, this modelling of the system is a clever way to avoid 
the fermionic sign problem. But it is more than a theo- 
retical trick in the case of ultra-cold atoms since it can be 
achieved experimentally by trapping atoms at the nodes 
of an optical lattice in the tight-binding regime [2^ . The 
Bethe-Peierls zero range model is another commonly used 
way of modelling the unitary regime: pairwise interac- 
tions are replaced by contact conditions imposed on the 
many-body wave function [H, [H, M, M M, 21,21,2^- 
This model is very well adapted to analytical calculations 
in few-body problems [13, [111 but can also be useful to 
predict many-body properties like time-dependent scal- 
ing solution [U, the link between short range scaling 
properties and the energy of the trapped gas [3l|, and 
hidden symmetry properties (35j of the trapped gas. 

However, there is to our knowledge no general rigorous 
result concerning the equivalence between the discrete 
(Hubbard model) and the continuous (Bethe-Peierls) 
models for the unitary gas. As a crucial example, one 
may wonder if there is any few- or many-body bound 
state in a discrete model at the infinite scattering length 
limit. This is a non-trivial question, since the infinite 
scattering length corresponds to an attractive on-site in- 
teraction in the discrete model. 

In this paper, we address this question for two and 
three fermions in a cubic box with periodic boundary 
conditions, when the interaction range tends to zero for 
a fixed infinite value of the scattering length. Our results 
for the equivalence of the lattice model and the Bethe- 
Peierls approach are analytical for two fermions but still 
rely on a numerical step for three fermions. In this few 
body problem, the grid spacing can however be made 
very small in comparison to the grids currently used in 
Quantum Monte Carlo many-body calculations, thus al- 



2 



lowing a more precise study of the zero lattice step limit 
and a test of the linear scaling of thermodynamic quanti- 
ties with the grid spacing used in . Our computations 
also exemplify the remarkable property that short range 
physics of the binary interaction does not play any signif- 
icant role in the unitary two-component Fermi gas, and 
the fact that the Bethe-Peierls model is well behaved for 
three equal mass fermions. 

Our model is the lattice model used in the Quantum 
Monte Carlo simulations of |24| . It has already been de- 
scribed in details in Refs. [3a, [33| so that we recall here 
only its main features. The positions of each particle i 
are discretized on a cubic lattice of period b. The Hamil- 
tonian contains the kinetic term of each particle, p^/2m, 
such that the plane wave of wave vector k has an energy 



2m 



(1) 



Here the wave vector is restricted to the first Brillouin 
zone of the lattice : 



keV= [-7r/6,7r/6p 



(2) 



We enclose the system in a cubic box of size L 
with periodic boundary conditions, so that the com- 
ponents {ka}ae{x,y.z} of k are integer multiples of 
2-k/L. In what follows we shall for convenience re- 
strict our computations to the case where the ratio 
L/h = 2N -f- 1 is an odd integer, so that ka = 2Trna/L 
with Ua e {~N,~N + 1, . . . ,N}. The Hamiltonian also 
contains the interaction potential between opposite spin 
fermions i and j, which is a discrete delta on the lattice: 



go 



(3) 



In [37j the matrix elements of the two-body T-matrix 
(k|r(i? -t- iO+jjk') for an infinite box size are shown to 
depend only on the energy E, not on the plane wave 
momenta, which would imply in a continuous space a 
pure s-wave scattering. The bare coupling constant t/o is 
then adjusted in order to reproduce in the zero energy 
limit the desired value of the s-wave scatte ring length a 
between two opposite spin particles [ssl. [36l. [3^: 
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where 



12 /"^/^ 

K ^ — I de\n{l + 1/cos^ 61) = 2.442749 ... (5) 

Jo 

may be expressed in terms of the dilog function, and 
g = Airh^a/m is the usual effective s-wave coupling con- 
stant. From the calculated energy dependence of the T- 
matrix, one may also extract the effective range of the 
interaction in the lattice model; r,, is found to be pro- 
portional to the lattice period, rg ~ 0.337 6 [H], and the 
limit & ^ is equivalent to the limit of both zero range 



and zero effective range for the interaction [331 ■ men- 
tioned in the introduction, this is the desired situation to 
reach the unitary limit when \a\ — oo. 

We first solve the problem for two opposite spin 
fermions in the box, in the singlet spin state |s) = (| ti 
) ^ I iT))/v^> by looking for eigenstates of eigenenergy 
E with a ket of the form |s) (g) 1^. We restrict to the 
case of a zero total momentum [39[, so that the orbital 
part 10) may be expanded on |k, — k) = |1 : k) (g) |2 : — k), 
where 1 1 : k) is the normalized ket representing particle 1 
with wave vector k. The corresponding wavcfunction is 
(r|k) = e^^'^ / L^^"^ . Schrodinger's equation then reduces 
to: 



(2ek-£;)(k,-k|, 
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(6) 



where the last term does not depend on a common po- 
sition r of the two particles. A first type of eigenstates 
corresponds to (r, r|(/)) = 0: these eigenstates have a zero 
probability to have two particles at the same point, and 
are also eigenstates of the non-interacting case. An exam- 
ple of such a state with the correct exchange symmetry 
is given by the wavcfunction: 



(ri, r2) oc cos 
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(7) 



We are interested here in the states of the second type, 
what we call 'interacting' states, such that (r,r|0) 7^ 0. 
Treating the interacting term in Eq. ([5]) as a source term, 
one expresses |(/)) in terms of (r, r|</)) and a sum over k. 
Projecting the resulting expression onto |r,r) leads to a 
closed equation (now E ^ 2ek): 
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The resulting implicit equation for E, of the form u(E) = 
0, where u{E) is the left hand side of Eq.®, is then read- 
ily solved numerically; to this end, one notes that u{E) 
has poles in each E = 2ek, and that it varies monotoni- 
cally from —00 to +00 between two poles, so that u{E) 
vanishes once and only once between two successive val- 
ues of 2ek. In Fig[l] we show for \a\ = 00 the first low ly- 
ing eigenenergies as functions of the lattice spacing; one 
observes a convergence to finite values in the b/L 
limit, with a first correction scaling as b/L. A rewriting 
of the implicit equation for E that will reveal convenient 
in the b = limit is: 



ttL 
a 



{2Tihf 



E ^ \E-2ek 2ek 



where the function C{b) is defined by: 



C(6) = 



{2TihfL 
2m 



d^k 



1 1^1 



L3 ek 



C'ib), 
(9) 

, (10) 
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FIG. 1; First three eigenenergies for the interacting states 
of two fermions in a box of size L for an infinite scattering 
length in the lattice model, as functions of the lattice period b. 
The total momentum of the eigenstates is fixed to zero. The 
computed eigenenergies are given by the plotting symbols, in 
units of Eo — (27r/i)^/2mL^; the straight lines are linear fits 
performed on the data with h/L < 2 x 10~^. 



where /lo-(q) = [exp(— g^a^) — and the integration 

domain is X = [-1/2, l/2]3. The desired identity C(0) = 
Cbp results from the fact that one can exchange the ct = 
limit and the summation over n in the above equation 

BUI. 

In the lattice model, it is possible to show analytically 
that the spectrum of the two-body problem for an infinite 
scattering length is bounded from below in the 6 — > 
limit. Since go < for |a| = cx), there exists at least one 
non-positive energy solution, by a variational argument. 
One then notes that the right hand side in Eq.® is a 
strictly decreasing function of E over ] — oo, 0[, that tends 
to — oo in i? = 0~, so that at most one negative energy 
solution may exist. Furthermore one can show that the 
& — !■ limit of the right hand side tends to +oo when 
E ^ ^oo 144] . whence this negative energy solution is 
finite m. 

We now turn to the problem of three interacting 
fermions in the box. Schrodingcr's equation is obtained 
without loss of generality by considering the particular 
spin component (1: t ; 2: t ; 3: J,), so that the interaction 
takes place only among the pairs (1, 3) and (2, 3), and in 
the lattice model one obtains: 



and has a finite limit for 6^0 which is given by C(0) ~ 
8.91364. 

We now briefly check that the 6 = limit in Eq.® 
coincides with the prediction of the Bethc-Pcicrls model, 
which is a continuous space model where one replaces the 
interaction potential by the following contact conditions 
on the wavefunction [ll, [13, [H, [H, [sO, [Ml, [H, i!] : there 
exists a function S'(R) such that, 



(ri,r2) = ^(R) 



1 



0(r), 



(11) 



where r = |ri — — > is the distance between the two 
particles and the center of mass position R = (ri -I- r2)/2 
is fixed. At positions ri 7^ r2, the wavefunction solves the 
free Schrodinger equation. Using this model we arrive at 
an implicit equation for the energy of an interacting state 
exactly of the form obtained by taking the 6 = limit in 
Eq.®, except that the constant C(0) in the right hand 
side is replaced by [4Ci| : 
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'0(ri,r2,r3) = 0. 



(14) 

We restrict to a zero total momentum modulo 27r/6 along 
each direction [39| : using the fcrmionic antisymmetry 
condition for the transposition of particles 1 and 2, we 
express the part of Eq. p4p involving the interaction in 
terms of a function of the position of a single particle: 



?/'(ri,r2,ri) = f{v2 - ri) 
V'(ri,r2,r2) = -/(ri - v^). 



(15) 
(16) 



We then project Eq. p^ on the plane waves in the box, 
which leads to: 



(ki,k2,k3|V') = 



„ ;:mod 

t/0 "ki+k2+k3,0 



(/k.-/kj (17) 



where (5™°^ is a discrete delta modulo 27r/6 along each 
direction, and where the Fourier transform of /(r) is de- 
fined as: 



Cbp = hm 

<T->0 



E 



(12) 



We expect the identity Cbp = C(0) from the general re- 
sult that the Bethe-Peierls model for the two-body prob- 
lem reproduces the zero range limit of a true interaction 
potential [H, [4l| . It is however instructive to check this 
property explicitly for the lattice model. One can show 
that: 

Cbp - C(0) = lini ^ j d^u [h„{n + u) - h,{n)] 

(13) 



/k=(k|/) = -^ E exp(-zk.r)/(r). (18) 



re[0,-L[- 



Replacing /(r) in the right-hand side of this equation 
by its expression in terms of (ki, k2, 'k.al'fp) deduced from 
Eq. p^ . we obtain a closed equation for /k: 
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/k = /k ^ flk.q - ^ Ok.q/q 



qeD qeP 

where we have introduced the matrix: 

1 



flk,. 
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ffkH 



(19) 



(20) 
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and for an arbitrary wavevector u, [u]fbz denotes the 
vector in the first Brillouin zone that differs fi'om u by 
integer multiples of 27r/& along each direction. The eigen- 
values E of the three-body problem are such that the 
linear system (jl9p admits a non-identically vanishing so- 
lution /k, that is the determinant of this linear system is 
zero. Note that from Eq. fTO)) . one has /(O) oc X]qe_D /k = 
0, a consequence of Pauli exclusion principle. 

For \a\ = GO, we have computed numerically the first 
eigenenergies of the system, by calculating the determi- 
nant as a function of E. In Fig[2]we give these eigenen- 
ergies as functions of the ratio b/L. A rapid convergence 
in the zero-6 limit is observed, with a linear dependence 
in b/L. 

This rapid convergence illustrates the fact that equal 
mass fcrmions easily exhibit universal properties, as re- 
vealed by experiments; here b plays the role of the fi- 
nite Van der Waals range of the true potential [given by 
{mCe/h^y/'^, where Cq is the Van der Waals coefficient], 
and L is of the order of the mean interparticle distance 
in a real gas. As an example, for ^Li atoms 6 ~ 3 nm and 
in experiments for the broad Feshbach resonance in the 
s-wave channel at ^ 830 G the atomic density is of the 
order of 10^"^ cm~^, so that the ratio b/L is of the order 
of 10~^ which is well within the zero-& limit. 

The absence of negative three-body eigenenergies in 
the unitary limit can be obtained numerically very effi- 
ciently through a formal analogy between Ea. p9)) and a 
set of rate equations on fictitious occupation numbers of 
the single particle modes in the box. Assuming E < 0, we 
note Ilk tfic fictitious occupation number in the mode k 
and Fk^q = gociq^k/L^ the transition rate from the mode 
k to the mode q. From Eg. ([^0)1 . one obtains the property 
Fk_,q = Fq^k, and the rate equation can be written as: 



dllk 

dt 




(21) 



The symmetric matrix M{E), which defines the first or- 
der finear system in dU/dt = M{E)Il, has the 
following properties: 1) its eigenvalues are non-positive, 
since it is a set of rate equations; 2) its eigenvalues are 
decreasing function of the energy E^ which can be de- 
duced from the fact that dM{E)/dE is a matrix of rate 
equations and obeys property 1), and from the Hellman- 
Feynman theorem; and 3) eigenmodes of Ea. ipTI) with 
an eigenvalue equal to —1 correspond to solutions /k of 
Eg. p^ with Ilk = /k cxp(— t). Therefore, in order to 
check that there is no non-zero solution of Eq. (fT9|) for 
iJ < 0, it is sufficient to check that all eigenvalues of 
M(E = 0) arc strictly larger than —1. 

We have computed the lowest eigenvalue toq of the 
matrix M{E = 0) as a function of the ratio b/L. A fit 
of mo as a function of b/L suggests limt,_to "^o — ~1- 
To better see what happens in the zero b/L limit, we 
note that having vtlq > — 1 is equivalent to having (mo 4- 
l)/5o < Oj or more simply (mo -f l)/{b/L) > 0. We have 
thus plotted in Fig[3]the ratio (mo + l)/{b/L), which is 
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FIG. 2: First eigenenergies of three fermions in a box of size 
L for an infinite scattering length in the lattice model, for 
a zero total momentum. The computed eigenenergies (dia- 
monds) are given in units of Eq — /2mL^ for different 
values of the lattice period h. For functions /(r) invariant by 
reflection along x, y, z and by arbitrary permutation of x, y, z 
we have computed the eigenenergies down to smaller values 
of h/L. The straight lines are a linear fit performed on the 
data over the range h/L < 1/15, except for the energy branch 
E ~ 2.89-Eo which becomes more slowly linear than the other 
branches. The eigenenergies predicted by the Bethe-Peierls 
model are given by stars in 6 = 0. 



seen to tend to a positive value for 6 — *■ 0, ~ 1.085, with 
a negative slope; this excludes the existence of negative 
eigenenergies for the three fermions at infinite scattering 
length even in the small b limit [i^ . 

In a last step, we compare the results of the lattice 
model to the predictions of the Bethe-Peierls approach 
for three fermions in a continuous space, which was shown 
to be a successful model in free space [l^,!!^] and in a har- 
monic trap at the unitary limit |32| . For this purpose, we 
introduce the function which is the Fourier transform 
of the regular part of the wave function as |ri — rsj ^ 0: 



F(R) = lim 



ri/) (R+ -,0,R- - 



(22) 



where we have used the translational invariance. By re- 
producing a calculation procedure analogous to what we 
have done for the lattice model, we obtain the following 
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FIG. 3: Quantity (mo + l)/(6/L) as a function of the lattice 
period b. Here mo is the lowest eigenvalue of the matrix M{E) 
defining the linear system Ea. (|21|> . for E = and for an 
infinite scattering length. The fact that mo + 1 > shows that 
there is no negative eigenenergy for the three fermions, see 
text. The symbols are obtained from a numerical calculation 
of mo. The solid line is a linear fit over the range h/L < 
1/29, not including the point with h/L — 1/81: for this point, 
the matrix M has more than half a million lines so that mo 
was obtained by a computer memory-saving iterative method 
rather than by a direct diagonalisation. 



infinite dimension linear system: 



L3 

— i^k = i^k 

9 



where the wavevectors k and q now run over the whole 
space (27r/L)Z'^, and: 



A 



1 



k,q 



E — e]i — Eq — Ck+q 



(24) 



The similarity between the structure of (|19p and ([23|) is 
apparent. Numerically, at \a\ = oo, we have verified the 
convergence between the two models as 6 in Eg. (fT9|l . 
see Figl^l Analytically, one can even formally check the 
equivalence between the two sets of equations (|19p and 
([23|) : First, we eliminate the integral of 1/ek between 
(|4|) and (fTO|) . to express 1/go in terms of l/g and C{b). 
Second, we replace 1/go by the resulting expression in 
Ea. p9|) . Third, we take the limit 6-^0: we exactly re- 
cover the system ([23|) . Hence, if the eigenenergy E and 
the corresponding function / in the lattice model have 
a well defined limit for 6 = 0, this shows that the limit 
is given by the Bethe-Peierls model. Of course, the real 
mathematical difficulty is to show the existence of the 
limit, in particular for all eigenenergies. This property 
is not granted: For example, the present lattice model 
generalized to the case of a | particle of a mass dif- 
ferent from the mass m of the two t particles leads, for 
a large enough mass ratio m/m^, to a three-body energy 
spectrum not bounded from below in the 6 = limit, 
even though the Pauli exclusion principle prevents from 
having the three particles on the same lattice site [47| . 

In conclusion, we have computed numerically the low 
lying eigenenergies of three spin- 1/2 fermions in a box, 
interacting with an infinite scattering length in a lattice 
model, for a zero total momentum and for decreasing val- 
ues of the lattice period. Our results show numerically 
the equivalence between this model and the Bethe-Peierls 
approach in the limit of zero lattice period. This is re- 
lated to the fact that the eigenenergies E are bounded 
from below in the zero lattice period limit 6 — > 0, more 
precisely E > 0. Such a convergence of the eigenstates of 
fermions in a lattice model towards universal states when 
6 ^ is a key property used in Monte Carlo simulations 
at the iV-body level (al, [13, [H] . 

We thank F. Werner for interesting discussions on 
the subject. Laboratoire de Physique Theorique de la 
Maticre Condensce is the Unite Mixte de Recherche 7600 
of Centre National de la Recherche Scientifique (CNRS). 
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